A simple model for calculating magnetic nanowire domain wall fringing fields 
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We present a new approach to calculating magnetic fringing fields from head-to-head type domain walls in 
planar magnetic nanowires. In contrast to calculations based on micromagnetically simulated structures the 
descriptions of the fields are for the most part analytic and thus significantly less time and resource intensive. 
We begin with an intuitive picture of domain walls, which is built upon in a phenomcnological manner. The 
resulting models require no a priori knowledge of the magnetization structure, and facilitate calculation of 
fringing fields without any free parameters. Comparisons with fields calculated using micromagnetic methods 
show good quantitative agreement. 



I. INTRODUCTION 

Magnetic domain walls are boundaries between ar- 
eas of differing magnetization direction, with magneti- 
zation rotating through the domain wall width. In pat- 
terned magnetic nanowires, domain walls separate op- 
posite magnetizations, either in a head-to-head (magne- 
tizations point towards the wall) or tail-to-tail (magne- 
tizations point away from the wall) configuration. The 
domain walls form in two characteristic shapes - trans- 
verse or vorte^^l. The former is more common in wires 
with smaller cross sections. The magnetization struc- 
tures of typical domain walls are shown in Figure [TJ as 
calculated via micromagnetic methods. The converging 
or diverging magnetization of a domain wall results in an 
associated effective magnetic monopole moment. Thus, 
strong magnetic fields are found in close proximity to the 
domain walls. Directly above the walls these fields are 
directed out of the nanowire's plane. 

The characteristics and dynamics of nanomagnetic do- 
main walls have been the subjects of much research, 
e.gP and are key topics in the burgeoning field of 
spintroniciP^. The fields of atom optics and atom chips 
have also experienced rapid expansion and continue to 
strive towards improved techniques for producing arenas 
for atomic physics experiments and quantum information 
processing^!. Within these research fields nanomagnetic 
domain walls have a number of possible applications, 
such as an atom mirrop^ or mobile atom traps^H A fac- 
tor common to all these areas of study is the need for an 
accurate knowledge of the magnetic fringing fields gener- 
ated by domain walls. In particular, for atomic physics 
applications a knowledge of the fields for distances of 
around 100 nm or greater is desired; at distances closer 
than this the attractive van der Waals interaction begins 
to modify significantly the atomic potential. 

Rigorous micromagnetic methods currently allow ac- 
curate computation of magnetization configuration, and 
hence fields, at all distances. However they are resource 
and time intensive and require detailed understanding of 
micromagnetics. We present a more direct method of 
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FIG. 1. (a) Magnetization structures of typical head-to-head 
domain walls as calculated via micromagnetic methods us- 
ing OOMMF^. The top image shows a transverse-type wall 
(200 nm x 5 nm wire cross section), the bottom image is 
a vortex-type wall (200 nm x 15 nm wire cross section). 
Shading represents the longitudinal component of the nor- 
malized magnetization, m. Arrows represent the in-plane 
magnetization direction, (b) Divergence of the magnetiza- 
tion. Dark (light) areas signify positive (negative) V.M, or 
effective 'north' ('south') poles. Cross sections are as in (a). 

field calculation - the 'monopole model', which affords 
considerable perspicacity and utility, particularly in the 
regimes of interest for applications of domain walls in 
atomic physics experiments. In Section [TT] we outline the 
models used, extending to greater complexity in a phe- 
nomenological manner. We discuss how emulating the 
detailed magnetization structure of the wall is not a good 
way of quickly and accurately computing these fields; the 
models we present do not demand this and can be used 
without free parameters. Analysis of the models' accu- 
racy is described in Section [ITT] 
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II. MODELS 

Illustrations of the four different models considered are 
shown in Figure [2] The basic 'monopole model' is de- 
scribed in Section lll Al and more detailed models are de- 
veloped in subsequent sections. Derivations are provided 
in Appendix |A"[ 



where f is the unit vector associated with r. Note that 
B and H differ only by a factor of /Uo- 

As will be shown in Section [mi this model is an excel- 
lent approximation at distances that are large compared 
to the width of the nanowire. We now extend the model 
to emulate more faithfully the shape of a real domain 
wall. This will significantly extend the region over which 
the model accurately reproduces the shape of real mag- 
netic fields. 



A. Monopole Model 



Although Maxwell's equations preclude the existence 
of magnetic monopoles . they have been posited and ob- 
served as quasiparticles^ES. We will use them as a the- 
oretical construct; at large distances, nanomagnetic do- 
main walls can be well approximated as point magnetic 
charges. This point acts as a source or sink of magnetic 
field lines - a monopole. This prescription of an effective 
charge is a treatment followed elsewhere, see e.g! 15 * 17 * ^ . 
The charge is then given by 



B. ID Domain Wall 

The first extension to the model is to consider the wall 
as a ID object, i.e. a line of charge across the width of 
the wire as illustrated in Figure |2b| The corresponding 
expressions for the field are given in Eqs. @-([5]). These 
expressions will be shown to offer a significant improve- 
ment in the accuracy of field calculations compared to 
the simple monopole model. 



C. 2D Domain Wall 



q m = 2[i M s wt, 



(1) 



where M s is the saturation magnetisation, /zo is the per- 
meability, w is the wire width and t is the wire thickness. 
This expression has previously been derivecP- the details 
are provided in Appendix |Aj The magnetic flux density 
('magnetic field' from hereon) for a charge q m has an in- 
verse square dependence on distance: 



B(f) = 
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We now extend our model to consider the domain wall 
to be 2D as a more exact representation of it as an ex- 
tended object. Modelling a domain wall with finite size 
in z does not give significant improvement; the nanowire 
thickness is small compared to its width. From hereon 
we consider domain walls with zero size in z, but finite 
extent in x and y, as illustrated in Figure [2c] The corre- 
sponding expressions for the magnetic field are shown in 
Eqs. 

Note that Equations ^ and Q are identical under 
an exchange of {x, w} with {y, sf, as expected due to 
symmetry. 
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D. Triangular Domain Walls 



shape associated with the two wall types; whilst the poles 



Considering the magnetic pole distributions, V.M, 
shown in Figure ([!]), there is a clear difference in the 
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(a) A domain wall represented as a region of 
magnetization reversal. In the limit s — > the 
wall is a 2D object of area wt. 



(b) A ID line of magnetic charge. Here the wire 
length is into the page. 





(c) A 2D area of magnetic charge with finite size 
along the wire length. 



(d) A triangle of magnetic charge. This shape is 
applicable for transverse walls. 



FIG. 2. Illustrations of the geometries of the four different models of the magnetic charge associated with a domain wall. 



of a vortex wall are contained within an approximately 
rectangular section of the nanowire, the pole distribution 
of a transverse wall is distinctly triangular. The final ex- 
tension of our model is to incorporate this characteristic 
shape, so as to better emulate the field from a transverse 
wall. This situation is illustrated in Figure |2d| 



Equations (|A6|) and ( A9 ) can be modified to incorpo- 



rate a triangular shape. Unfortunately it was not possible 
to derive closed expressions for B z or \B\ as with previous 
models. Whilst numerical integration is quicker than us- 
ing micromagnetically simulated structures for the same 
level of resolution, it is significantly slower than previous 
analytic expressions. Instead, the computation can be 
sped up considerably by utilizing the result for a rectan- 
gular wall. One can approximate the shape of the trian- 
gular wall by a series of rectangles and sum the result- 
ing field from each. Computing the field over 10 6 points 
takes around 2 seconds per rectangle. Negligible loss of 
accuracy was observed using a triangle divided into 40 
rectangles. 

Whilst the models presented thus far emulate the basic 
shape of the domain walls' volume charge distributions, 
they do not take into account the positive and negative 
regions of edge charge present in both transverse and 
vortex walls (Figure . As can be seen in Figure [3| at 
extremely short range the fields produced by the domain 
walls mimic this complex charge structure. A number of 
efforts were made to emulate these more subtle features; 



trapezoidal walls, dominating edge charge regions and 
simple spatial variations in charge density /polarity were 
all investigated. However, none of these approaches pro- 
duced overall improvements in accuracy and therefore we 
do not describe them in detail here. Thus we conclude, 
in agreement with previous findings^, that the overall 
shape of a domain wall's volume charge is by far the 
most important feature at all but the very shortest dis- 
tances. Whilst the accuracy of the models in the near 
field could undoubtedly be improved by incorporating a 
detailed magnetization structure, it is precisely this we 
are trying to circumvent with the development of these 
simple models. 

In the 2D models we use the approximation s = w. It is 
possible to optimize the value of s through a comparison 
with micromagnetic simulations, however this approach 
is contrary to the aim of these models and significantly 
detracts from their utility. We will see later that us- 
ing an optimized value of s confers only a very small 
improvement in accuracy. Unless explicitly stated oth- 
erwise, s = w will be used throughout. We note for a 
triangular wall (as with a rectangular one) that the dis- 
tribution of charge across % is independent of s: 

dq m (x^) = 4M s /ini/s dA — 4M s fJ,QtxN /w dx^- (9) 

This is equivalent to the fact that all triangles have cen- 
ters of mass at barycentric coordinates (1/3, 1/3, 1/3). 
The field maximum according to our analytic model is 
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thus found at (x, y) = {w/3, 0) when sufficiently far from 
the wall. Note also integrating d<7 m (xN) with respect to 
x N yields Equation (IA4I). 



III. COMPARISON OF MODELS 

To assess the accuracy of the models presented com- 
parison was made with results from calculations based 
on micromagnetic simulations. We use a proprietary mi- 
cromagnetic code that solves the Landau-Lifshitz-Gilbert 
equation of motion and quasistatic Maxwell equations 
within a finite element/boundary element frameworlP^. 
We simulate 8.4 urn long nanowire structures, discretized 
within the domain wall region into tetrahedral meshes 
with a 5 nm characteristic size. Physically appropri- 
ate domain wall structures are introduced by imposing 
simple bi-domain states and then allowing relaxation to 
equilibrium. The magnetic field profiles above the do- 
main walls are calculated analytically from the e quiva- 
lent dipole charges on the nanowires ' surfacefPE] Mag- 
netic fields created by effective magnetic charges at the 
nanowire ends are subtracted from the data using the 
simple point monopole approximation, which, as we will 
show later, is extremely accurate in the far-field. Fig- 
ure [3] shows the increase in compexity of the field shape 
at short distances. Because of this complex shape the 
models presented here will always experience inaccuracy 
when calculating fields at very short distances. 

We consider the fields calculated by all four models, for 
a transverse-type domain wall in a wire with w = 200 nm 
and t = 5 nm. A saturation magnetization, M s , of 
8.6 x 10 5 A/m was used throughout^. The models were 
examined over a 1 |xm x 1 urn x 1 |xm cube divided 
evenly into 10 6 points, with the domain wall centered at 
the bottom of the cube. This is representative of dis- 
tances within which atomic physics applications of do- 
main walls aim to work; outside this region the field is 
less than 1 G. 

The following analysis examines only the field magni- 
tude for the sake of brevity. Very good accuracy was also 
observed for the field direction and is discussed briefly to- 
wards the end of the section. A summary of the different 
figures of merit is provided in Appendix [5] Initial com- 
parison was made by considering the maximum field at 
a given height, shown in Figure [4] This is an important 
quantity in atomic physics applications of domain walls, 
and shows well how the models scale with distance. 

From Figure [4] it can be seen that at heights greater 
than 50 nm there is very good agreement (< 10% error) 
with micromagnetic simulations from all but the simplest 
'monopole' model. There is a stark improvement moving 
to the ID case - for z < 300 nm there is an order of mag- 
nitude improvement in the error. Above this height all 
the models converge with the micromagnetic simulations, 
having an error of less than 5%. Extension beyond the 
ID model does not afford significantly more accuracy. 

A more thorough analysis of the models examines the 



field over the entire region of space, so the shape of the 
field is tested. Figure [5] shows the RMS error over all 
points at a given height for the different models consid- 
ered. We see that there is a big improvement in adopting 
a triangular shape; more accurately imitating the shape 
of the domain wall produces a more accurate field shape. 
Typical values of the RMS and mean percentage errors 
are also provided in Appendix [B] 

As can be seen in Figure [5j the accuracy of all the 
models is much worse at small distances, due to the dra- 
matic change in the shape of the fields. An examina- 
tion of e.g. Figure [3] suggests highly localised regions of 
magnetic charge for both wall types. However, as pre- 
viously mentioned, improved accuracy is not gained by 
using a simple representation of these charge distribu- 
tions. Whilst very good accuracy is achieved for all wire 
geometries at heights >100 nm, some spintronics applica- 
tions require knowledge of the field very close to the wire, 
e.g. read/ write heads can be located at flying heights of 
~10 nrrP^. This is within the regime where these analytic 
models break down. 

These initial figures of merit are strong indicators of 
the utility of the model we present, allowing non-expert 
users to calculate fields at all but the closest distances 
with a very good level of accuracy. Whilst mean errors 
indicate well the overall accuracy of models, there may 
be small regions of space where very large errors exist. 
As such, maximum errors are also provided: Em&x is the 
maximum % error in the field, is the maximum % 

error above 100 nm. Ai?Max and A5^° are the corre- 
sponding absolute errors in the field. A summary of the 
figures of merit is given in Appendix [B] 

The best overall representation of the accuracy of the 
models is achieved by examining the distribution of the 
error. This is shown for the 200 nm x 5 nm wire, using 
the triangle model, in Figure [6] At all heights the large 
majority of points have an error < 20%. The error di- 
rectly above the triangle barycenter (white line) is quite 
high within the points at a given height; this is the point 
closest to the wall for a given height. 

In an effort to present a figure of merit independent 
of sample size and distance from the wall center, we now 
perform analysis over regions of specific field strength. 
This removes biases due to changes in the volume of space 
examined, and will also be independent of nanowire size, 
whilst still assessing the fidelity of the models to the true 
field shape. The regions analyzed are shells centered on 
field isosurfaces. Whilst for some applications the field at 
a given distance is the focus, there are a number of appli- 
cations, e.g. atom trappingES, where the working regime 
is defined by the field strength itself. The following anal- 
ysis illustrates the error across these different regimes of 
field. We compute the error over a region where, accord- 
ing to micromagnetic simulations, 0.9i?o < \B\ < 1.1-Bq, 
for some Bo . The results of this analysis are shown in Fig- 
ure [7J A familiar trend is observed; the error increases for 
larger fields (shorter distances), and the error decreases 
for models which more faithfully replicate the shape of 
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Transverse Vortex 




FIG. 3. Fields calculated via micromagnetic simulations, (a) and (b) are at z = 200 nm, whilst (c) and (d) are at z — 12.5 nm; 
z = is located in the center of the wire's thickness, (a) and (c) correspond to a transverse-type domain wall (cross section = 
200 nm x 5 nm), whilst (b) and (d) correspond to a vortex-type wall (cross section = 200 nm x 15 nm). At short distances the 
field adopts a complex shape, indicative of the magnetic charge distribution. Further away the field is significantly smoother 
and smaller in magnitude. Note in the far field the shape is skewed in x for the transverse domain wall. 




FIG. 4. Maximum field at a given height above a transverse 
domain wall (200 nm x 5 nm wire) as calculated by the var- 
ious models. Note the good agreement at >50 nm for all 
models except the simple monopole model. 



FIG. 5. RMS accuracy of the field over all points at a given 
height above a transverse domain wall (200 nm x 5 nm wire). 
Note the monopole model is very inaccurate at shorter dis- 
tances, and that incoporating a triangular shape confers sig- 
nificantly increased accuracy. 



real domain walls. For very large fields the differences be- 
tween models become less clear. This is expected as the 
shape is not accurately reproduced by any of the models 
in the near field. Again there is convergence at low fields 
as the approximation of the domain wall as a point object 



becomes increasingly valid. The biggest benefit of using 
a more complex model is seen in the range ^10-100 G. It 
is precisely this region where the shape of the field begins 
to reflect the asymmetry of the domain wall shape, but 



G 



Percentage 
if Points 



10 



10" 



10 1 



Percentage Error 




FIG. 6. The distribution of the percentage error in the field 
at all heights as calculated by the triangular model for a 
200 nm x 5 nm nanowire. At each height the data are grouped 
into 100 evenly distributed bins. The white line shows the er- 
ror directly above the wall barycenter. 




Bo (G) 

FIG. 7. RMS accuracy of the field within regions cen- 
tered around field isosurfaces above a transverse domain wall 
(200 nm x 5 nm wire). For a given Bo, the error is averaged 
over points with field within 10% of Bq. 



does not show the complex behaviour observed at very 
short distances. 

Discussion up to this point has been based upon the 
field size. Similar analysis was also performed for the 
field direction. This has inherent difficulties as B x and 
By have zero points at all z. By examining the angle 
between the computed field directions we observe a good 
level of accuracy. E.g. consider the error in the field 
direction at a given height above the domain wall, using 
the triangular model. For a 200 nm x 5 nm wire the 
mean (maximum) error is less than 6° (20°) for all heights 
>80 nm. There is a significant reduction in accuracy at 
very short distance. This inaccuracy is, in some regions, 



stark. The reason for this is small regions of negative 
effective charge. For a transverse domain wall such a 
region is found at the point of the triangle. The resulting 
large error in B z is then responsible for the inaccuracy 
in field direction, which is prominent for z < 100 nm; 
higher than this there is excellent agreement, as with the 
field magnitude. 

The analysis presented here has been for one partic- 
ular wire geometry. To examine the flexibility of our 
approach to model nanowires with other dimensions we 
also considered five other wire geometries, detailed in Ta- 
bic [T] Note that the three geometries with larger wire 
thicknesses correspond to wires which host vortex-type 
domain walls. For these wires the triangular model is no 
longer appropriate since, as can be seen from Figure [T] 
the magnetic pole distribution for a vortex-type domain 
wall is rectangular in shape. A summary of the figures of 
merit for all models and wire geometries is provided in 
Appendix [Bj 



TABLE I. The six geometries of wire examined in comparisons 
with fields from micromagnetic methods. 



Label 


Width (nm) 


Thickness (nm) 


Wall type 


A 


100 


5 


Transverse 


B 


200 


5 


Transverse 


C 


400 


5 


Transverse 


D 


100 


15 


Vortex 


E 


200 


15 


Vortex 


F 


400 


1.5 


Vortex 



A general trend is observed that there is a loss of accu- 
racy with an increase in nanowire size, shown explicitly 
in Appendix [Bj This is due to larger wires having larger 
and more prominent 'near field' regions. Whilst the fig- 
ures presented only show one wire geometry, they are 
representative of all wire sizes; the distributions of error 
show a similar shape for all six wires considered. 

The analysis presented has used the assumption s = w. 
In Appendix |Bj we provide the RMS percentage error 
using an optimized value of s, -Erms- Comparing this 
with the RMS error using s — w, -Erms, it is clear that 
there is only a very small loss of accuracy when using 
this approximation. 



IV. CONCLUSIONS 

A series of analytic models of fields from domain 
walls in nanomagnetic structures have been derived and 
developed to incorporate the characteristic shape of 
transverse-type walls. The results obtained from these 
models were compared with micromagnetic simulations. 
Improved accuracy was observed for models incorporat- 
ing higher dimensionality in domain wall structure, re- 
flecting the fact that domain walls are not simple point 
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like objects. We have also shown that no a priori knowl- 
edge is required that might limit the usefulness of these 
models; assuming the wall width (along the wire) to be 
equal to the wire width results in negligible loss of accu- 
racy. 

Data from analytic models show better agreement at 
larger distances from the nanomagnetic structure - re- 
gions where the field is smaller. This is intuitively ex- 
pected as further from the domain wall the approxima- 
tion of it being a point object is increasingly accurate. At 
very short range more detailed structures are observed in 
micromagnetic simulations. The models presented in this 
paper do not reproduce this. 

Examples have been presented where an RMS error 
over an extended region of interest of less than 5% is 
achieved. For fields up to ^100 G an error of less than 
10% can be achieved. The maximum field for a given 
height is reproduced accurately by all but the simplest of 
models. Similar accuracy is also obtained when consid- 
ering the direction of the field. 

The models presented give an efficient and intuitive 
way of calculating the fields from domain walls in 
nanowires. The loss of accuracy compared to detailed 
micromagnetic simulations is small. The methods pro- 
vided are quicker and much more accessible than existing 
techniques. 
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Appendix A: Fringing field expressions 

Here we provide the derivations of the expressions for 
the magnetic field for all the presented models. The sym- 
bols used throughout are defined as follows: M s = satu- 
ration magnetization, /i = permeability of free space, w 
= wire/domain wall width, t = wire/domain wall thick- 
ness, s — domain wall length and {x, y, z} — spatial co- 
ordinates. 



The 'magnetic charge', 
is then given by 



associated with this volume 



9m = -Mo / V.M dV = -no / M.n dS. 



(A3) 



Here we employed the divergence theorem; ft is the unit 
normal of the surface element dS with S enclosing V. 

Assuming a discontinuous magnetization reversal the 
volume of charge tends to zero. This is illustrated in 
Figur e [2a| by considering s — > 0. Evaluation of Equa- 
tion (A3) yields the charge associated with a head-to- 



head wall, 



q m = 2fi M s wt. 



(A4) 



A tail-to-tail wall has a charge of — q m . Equation ( A4 ) has 



a form similar to other expressions for effective magnetic 
charge, c.f. e.gP^I. The magnetic field at a position r is 
then given, in direct analogy to Coulomb's law, by 



B(r) 



/ioM s wt . 
2n\rf 1 



(A5) 



2. ID Model 

We now extend the domain wall from a point charge 
to a line of charge. Infinitesimal elements of this charge, 
d<7 m , contained in a length dxN, are given by 



dq m = 2p M s t dxj, 



(A6) 



The contributions from across the entire wire width are 
summed in the following integral to give the components 
of magnetic field, 



Bi = 



where we define 



IM>M B t j w ' 2 r^ 
2?r J_ w/2 |rf 



dx> 



(A7) 



(r x ,r y ,r z ) = (x - x N , y - y N , z) 
the vector from an infinitesimal charge element located 
at (a;N,j/N,0) to the point under consideration (t/n = 
in this case). This then yields the expressions given in 
Eqs. (§-(§3 



1. Simple Monopole 



3. 2D Rectangle 



In the presence of a magnetic medium we have from 
Maxwell's equations: 



(Al) 



MoV. (H + M = 



where H is the magnetic field and M the magnetization. 
Considering the wall as an extended object (Figure 2a) 
with magnetic charge density p m we have in analogy to 
Gauss' law for electrostatics 



V.H = p m /no- 



(A2) 



The final extension to the model which we derive an- 
alytically is to consider the domain wall as a 2D object 



1 An alternative route to deriving these expressions is to use the 
magnetic scalar potential and integrating in an entirely analo- 
gous manner; the magnetic field is then simply the negative of 
the gradient of the scalar potential. Whilst entirely equivalent, 
within this paper we integrate expressions for the magnetic field 
directly throughout for the sake of clarity. 



in the xy plane, as in Figure 2c The wall size can now 
vary arbitrarily along the wire length, to have width s. 
We choose to fix the value of s to be w, the wire width. 
We rescale dg m to reflect the change in domain wall size. 
By analogy to Equation ( A6 ) we have 



dq m = 2fj, M s t/w dxNdjfr 



(A8) 



which is the infinitesimal element of charge associated 
with an area element dA —dx^dy^. Using this expres- 
sion the total field is then given by integrating over all 
dA. 



B, = 



w/2 rw/2 



HoM B trj 

w/2 J -w/2 2TTW |f| 3 



dyjsidxi 



(A9) 



Evaluation of these integrals was performed 
symbolically^, see also e.g.^. 



Appendix B: Figures of merit 

The following table shows all the figures of merit for the six different wire sizes according to the various models. 



TABLE II. A summary of the figures of merit. Erms is the RMS error over all points, E RMS is the RMS error with optimal 
s, .Em is the mean percentage error, Er^ms 3 is the RMS error in the maximum field for a given height, Em sx is the maximum 
observed percentage error in the field; E M *°° is this value limited to heights >100 nm, A£>Max is the maximum observed error 
in the field; AB^™ is this value limited to heights >100 nm. The labels A-F refer to the wire geometries detailed in Table [i] 
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